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Mode-I  crack  growth  under  conditions  of  generalized  plane 
stress  has  been  investigated.  It  has  been  assumed  that  near  the 
plane  of  the  crack  in  the  loading  zone,  the  simple  stress  components 
corresponding  to  a  central  fan  field  maintain  validity  up  to  the 
elastic-plastic  boundary.  By  the  use  of  expansions  of  the  particle 
velocities  in  the  coordinate  y,  and  by  matching  of  the  relevant 
stress  components  and  particle  velocities  to  the  dominant  terms  of 
appropriate  elastic  fields  at  the  elastic-plastic  boundary,  a 
complete  solution  has  been  obtained  for  e^in  t'ie  plane  °f  the 
crack.  The  solution  is  valid  from  the  propagating  crack  tip  up 
to  the  moving  elastic-plastic  boundary.  A  self-similar  crack 
nucleating  at  a  point  and  steady-state  crack  propagation  have  been 
considered  as  special  cases. 
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INTRODUCTION 


The  asymptotic  structure  of  quasi-static  near-tip  fields  of  stress 
and  deformation  for  a  growing  crack  in  an  elastic  perfectly-plastic 
material  has  been  discussed  by  several  authors.  A  detailed  discussion 
of  near-tip  fields  has  recently  been  given  by  Rice  (1982) ,  who  reviewed 
the  earlier  contributions  and  retrieved  them  as  special  cases  of  a 
general  formulation  for  materials  of  arbitrary  yield  condition  and 
associated  flow  rule.  In  general,  the  analytical  near-tip  results 
must  be  supplemented  by  numerical  calculations  to  obtain  certain 
arbitrary  functions  that  enter  in  the  asymptotic  considerations  of 
the  near-tip  fields. 

In  this  paper  analytic  expressions  have  been  obtained  for  the 
strain  in  the  plane  of  the  crack,  which  are  valid  from  the  propagating 
crack  tip  to  the  moving  elastic-plastic  boundary.  The  analytical 
approach  employs  expansions  of  the  particle  velocities  in  powers  of 
y  (the  distance  from  the  plane  of  the  crack) ,  and  yields  ordinary 
differential  equations  with  respect  to  x  for  the  coefficients.  The 
arbitrary  functions  that  enter  in  integrating  these  equations  have 
been  obtained  by  matching  the  fields  in  the  plastic  loading  zone  to 
the  dominant  terms  of  the  corresponding  elastic  fields  at  the  elastic- 
plastic  boundary. 

The  method  is  first  demonstrated  for  anti-plane  strain,  and  the 
results  of  Rice  (1968)  are  recovered.  Next,  the  case  of  generalized 
plane  stress  has  been  considered  in  detail.  It  has  been  assumed  that 
the  stress  components  for  a  centered  fan  field  (see  e.g.  Hutchinson, 
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1968) ,  which  satisfy  the  yield  condition  and  the  equilibrium 
equations,  are  valid  up  to  the  elastic-plastic  boundary  (at  least 
in  the  plane  of  the  crack) .  An  explicit  expression  has  been  obtained 
for  the  total  strain  e^Cx.O.t).  Self-similar  crack  growth  fields 
and  steady-state  fields  have  been  considered  as  special  cases.  The 
results  of  this  paper  are  particularly  suited  for  use  in  conjunction 
with  a  critical  strain  criterion. 
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2.  Governing  Equations 

The  x^-axis  of  a  stationary  coordinate  system  is  parallel  to 
the  crack  front,  and  points  in  the  direction  of  crack  growth.  The 
position  of  the  crack  tip  is  defined  by  x^  =  a(t).  A  moving  coordinate 
system,  x,y,z  is  centered  at  the  crack  tip,  with  its  axes  parallel  to 
the  x^»x2  an£*  x3  axes.  Relative  to  the  moving  coordinate  system  we 
also  define  polar  coordinates  r,0,  with  0=0  coinciding  with  the 
positive  x  direction.  The  geometry  is  shown  in  Fig.  1. 

Relative  to  the  stationary  coordinates  the  equilibrium  equations 


are 


Vu  " 0 

where  cr^  -  o^±  is  the  stress  tensor  (i,j,k  =  1,2,3),  =  a/3x±, 

and  the  summation  convention  applies. 

The  Huber-Mises  yield  criterion  may  be  stated  as 

i  Vij ' k! 

where  s..  is  the  deviator ic  stress  tensor  and  k  is  the  yield  stress 
ij 

in  pure  shear. 

The  rate  of  deformation  is  defined  by 


(2.1) 


(2.2) 


Eij  ‘  2<8jVi  +  8iV 


Here  is  the  (infinitesimal)  strain, 
time  rate  at  a  fixed  material  point,  v. 


the  superposed  dot  denotes 
is  the  velocity  field,  and 


(2.3) 
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Vi  *  Ui*  w^ere  ui  *s  c^e  displacement  field.  The  total  strain  rates 
are  defined  by 


Dij  “  °ij  +  °ij  * 


where  the  elastic  strain  rates  are  related  to  the  stress  rates  by 

"Ij  -  '•yV1  ■ 

and  the  plastic  strain  rates  are 


‘  A“u 


(2.4) 


(2.5) 


(2.6) 


In  (2.5)-(2.6),  E  and  v  are  Young's  modulus  and  Poisson's  ratio, 
respectively  and  A  is  a  positive  function  of  time  and  the  spatial 
coordinates. 

The  spatial  derivatives  with  respect  to  the  cartesian  coordinates 
in  the  stationary  and  moving  coordinate  systems  are  directly  related  by 


In  the  moving  coordinate  system  the  material  time  derivative  is 

(*)  -  8t  -  a3x  (2.8) 

where  a  -  da/dt  is  the  speed  of  the  crack  tip. 

3.  Anti-plane  Strain 

In  the  moving  coordinate  system  the  only  non-vanishing  displacement 

component  is  w(x,y,t).  The  non-vanishing  stress  components  are 

o _  ■  a  and  a  •  o  The  Huber-Mlses  yield  criterion  (2.2) 

xz  zx  yz  zy  • 


reduces  to 
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(a  )z  +  (a  )' 
xz  yz 

In  the  loading  region  ahead  of  the  crack  tip,  the  stresses  may 
be  written  as 


(3.1) 


a  m  k  cos0  ,  a  =  -  k  sinG  (3.2a,b) 

yz  xz 

It  may  be  verified  that  (3.2a,b)  not  only  satisfy  (3.1),  but  also  the 
equilibrium  equations  (2.1). 

For  small  values  of  0  (i.e.,  y/x  <<  1)  the  plastic  fields  in  the 
loading  zone  will  be  matched  to  the  dominant  terms  of  the  elastic 
fields.  In  polar  coordinates  R,i|>  centered  at  point  E  the  dominant 
terms  of  the  elastic  solution  on  the  elastic  side  of  the  elastic-plastic 
boundary  are  taken  in  the  general  form 


1  Km  sin|ip  ,  V  -  shear  modulus 


°Rz  "  Girtil*  K 

G.p/ 


iii  “‘"T1, 


Km  c°4" 

It  should  be  noted  that  the  center  of  the  elastic  field  is  not  taken 
to  coincide  with  the  crack  tip.  The  center  is  located  at  a  moving 
point  E,  whose  position  is  defined  by  x^  -  e(t).  The  geometry  is 
shown  in  Fig.  2. 

From  the  condition  that  the  elastic  field  should  just  reach  the 


(3.3) 

(3.4) 

(3.5) 


yield  condition  (3.1)  at  the  elastic-plastic  boundary,  we  obtain  by 
the  use  of  (3.4)-(3.5) 
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(3.6) 


where  R  *  R^  defines  the  radius  of  curvature  of  the  elastic-plastic 


boundary*  at  least  for  small  values  of 


Another  condition  is  that 


o„  (i.e.  the  shear  stress  component  in  the  R,ip  system)  should  be 
KZ 


continuous  at  the  elastic-plastic  boundary.  We  find  by  using  (3.2a,b): 


o  *  o  cosip  +  o  sinip  «  k  sin(ip-0)  (3.7) 

kz  y  z  xz 

Thus,  by  the  use  of  (3.4)  and  (3.6) 

sin|ip  =  sin(ip-0)  ,  i.e.,  0  =  •|ip  (3.8a,b) 

For  small  6  and  ip  this  result  implies  that  the  center  of  the  elastic 
field  E  is  located  at 


x  *  i  x  ,  and  thus  e(t)  -  a(t)  +  x  (t),  R  (t)  =  \  x  (t)  ,  (3.9a,b,c) 

b  2  p  2  p  p2p 

where  x^(t)  defines  the  x-coordinate  of  the  elastic-plastic  boundary  in 
the  plane  y  «  0.  Thus,  just  as  for  the  case  of  a  stationary  crack,  the 
elastic  field  is  centered  halfway  between  the  crack  tip  and  the  elastic- 
plastic  boundary. 

The  third  condition  requires  that  the  particle  velocity  is 
continuous  at  the  elastic-plastic  boundary.  In  the  moving  coordinates 
R»ip»  the  material  time  derivative  may  be  written  as 

(*)  -  at  -  e  (cosip  3R  -  3^)  (3.10) 

where  e  is  the  velocity  in  the  x^-direction  of  the  center  E  of  the 
elastic  field.  By  applying  the  operator  (3.10)  to  the  expression  for 
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w  given  by  (3.3),  we  obtain  at  R  *  R^ 

w  =  —  (R  +  e)  sin^  (3.11) 

P  W  P  2 

.  •  i  . 

where  the  relation  (3.6)  has  been  used.  Since  e  =  a  +  —  x 

2  P 

•  2.  • 

and  R  *  -r  x  ,  Eq.(3.11)  further  reduces  to 
P  2  p 

w  =  —  (a  +  x  )  sin0  (3.12) 

PUP 

where  (3.8b)  has  also  been  used. 

For  quasi-static  growth  of  cracks  in  anti-plane  strain,  the 

particle  velocity  in  the  loading  region,  has  been  derived  by  Rice  (1968) 

in  terms  of  the  position  of  the  elastic-plastic  boundary.  The 

corresponding  expression  for  the  total  strain  in  the  plane  of  the 

crack  in  between  the  crack  tip  and  the  elastic-plastic  boundary 

depends  only  on  a(t)  and  x  (t) .  This  interesting  result  suggests  that 

P 

it  should  be  possible  to  obtain  the  strain  independently  from  a 
closed  system  of  equations  valid  only  near  y  *=  0.  We  will  show  the 
derivation  of  the  relevant  equations  primarily  as  a  preliminary  to 
the  corresponding  analysis  for  the  more  complicated  plane  stress  case. 

As  point  of  departure  we  take  the  following  expansions  for 
y/x  «  1: 

w(x,y,t)  -  w^x.t)  y  +  w3(x,t)  y3  +  0(y5)  (3.13) 

A(x,y,t)  -  AQ(x,t)  +  A2(x,t)  y2  +  0(y1')  (3.14) 

The  strain  rates  follow  from  (3.13)  as 


1  s"l 


xz  ■  2  JT 1  *  °<y’> 


(3.15) 


D  *  —  w  +  —  w  v2  +  0(v**) 
yz  2  12  3  y  vy  1 

From  Eqs. (2.4)-(2.6)  we  find  by  using  (3.2a,b) 


(3.16) 


D  =  -  —  sin6cos6  -  Aksin0 

xz  2p  r 


(3.17) 


D  ■  -  —  sin20  +  Akcos0 

yz  2 p  r 


For  small  0  -  y/x  «  1,  equality  of  terms  of  order  unity  in  (3.16) 
and  (3.18)  yields 


Aok  •  2  »1 

Equality  of  terms  of  order  y  in  (3.15)  and  (3.17)  gives 

_  1  k  a  _  *  k  _  1  3wi 
2  p  x7  ox  2  3x 

By  combining  (3.19)  and  (3.20)  we  then  find 


(3.18) 


(3.19) 


(3.20) 


3wi  w,  ,  • 

_ 1  |  1  _  k.  a 

3x  x  “  ~  p  7 


(3.21) 


Integration  of  (3.21)  yields 


w  -  -  i  i  IM-)  +  ^ 
1  p  x  x  x 


(3.22) 


where  x  ■  x^(t)  is  the  length  of  the  plastic  zone  in  the  plane  of 
the  crack.  Evidently,  A(t)(y/x)  is  the  particle  velocity  at  the 
elastic  plastic  boundary  near  y»0.  For  small  values  of  0  it 
follows  from  (3.12)  that  A(t)  -  (k/p)(a  +  xp).  Integration  of  (3.22) 
over  tine  in  the  fixed  coordinates  yields  for  t  _>  t^: 


V 


ly'V0-0  +  t  rn  -  «n  +  £*)  dt  0.23) 

S  p 

where  is  the  time  that  the  elastic-plastic  boundary  arrives  at 

position  x,,  i.e.,  a(t  )  +  x  (t  )  =  x,  and  k/p  is  the  elastic  strain 
1  P  P  P  1 

at  the  elastic-plastic  boundary. 

Several  interesting  observations  based  on  (3.23)  have  been 

made  by  Rice  (1968).  As  discussed  by  Rice  a  criterion  of  critical 

plastic  strain  y  at  a  distance  p  ahead  of  the  crack  tip  (McClintock 
P  s 

and  Irwin,  1965)  may  be  employed  in  conjunction  with  (3.23)  to  compute 
the  required  x^(t)  for  quasistatic  extension  of  the  crack  over  a 
distance  a(t).  We  note  that  the  corresponding  value  of  Rp(t)  is 
given  by  (3.9c),  while  the  required  external  load  subsequently  follows 
from  (3.6). 
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4 .  Generalized  Plane  Stress 

For  generalized  plane  stress  the  Huber-Mises  yield  criterion  (2.2) 
becomes 

a2  +  a2  -  a  a  +  3a2  =  3k2  (4.1) 

x  y  x  y  xy 

In  the  loading  region  ahead  of  the  crack  tip  expressions  for  the 
stresses  which  satisfy  both  (4.1)  and  the  equilibrium  equations  are 

o  =  kcos30,  a  =  k(2cos30  +  3sin20  cos0) ,  a  =  -  ksin30  (4.2a,b, 
x  y  xy 

For  small  values  of  0  (i.e.,  y/x  <<  1)  the  field  in  the  plastic 
loading  zone  will  be  matched  at  the  elastic-plastic  boundary  to  the 
dominant  terms  of  a  corresponding  elastic  field.  For  the  elastic  field 
we  do,  however,  not  take  the  field  for  a  crack,  but  rather  that  for  a 
notch  with  -|p  as  radius  of  curvature  at  its  tip.  In  polar  coordinates 
R,i|/,  the  appropriate  Mode-I  stress  fields  are  given  by  Creager  and  Paris 
(1967)  as 

h 

°x  =  (inR  )  Ki{co4"  U-sin^  sin|^]  -  ^  cos|t/j)  (4.3) 

h 

°y  =  (2^)  Ki{coaf^  [1+sinf^  sirt|^]  +  cos|if/}  (4.4) 

°xy  =  (^r)  KI{sir 4*  COS¥  C08b  ~  2R  Si4}  (4*5) 

Note  that  the  tip  of  the  notch,  which  is  not  the  tip  of  the  crack 
nor  the  elastic-plastic  boundary,  is  a  distance  -|p  from  the  origin  E, 
as  shown  in  Fig.  3.  Just  as  for  the  Mode-Ill  case,  the  center  of  the 
elastic  field  E,  whose  position  is  defined  by  x^-  e(t)  ,  y^  ■  0,  is 
located  in  between  the  crack  tip  and  the  elastic-plastic  boundary 
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defined  by  x**x  (t).  For  generalized  plane  stress,  the  displacements 
P 

corresponding  to  (4.3)-(4.5)  are 

h 

u  =  Kjtcosl^  [k-1  +  2  sin2-|ijj]  +  £  cos|i{j}  (4.6) 

h 

v  =(h)  ly  Kx^ sin§^  t<+1  ~  2  cos2-^]  sin|^l  (4.7) 

where  k  =  (3-v)/(l+v). 

From  the  condition  that  the  elastic  field  should  just  reach  the 
yield  condition  at  the  elastic-plastic  boundary,  we  obtain  by  the  use  of  (4.1) 


(4.7) 


and  (4.3)-(4.5) 


[(sr)\]2  [1+3y-)2]  - 3k! 

p  p 


(4.8) 


where  R  =  R^  at  the  elastic-plastic  boundary,  at  least  for  small 

values  of  y.  Another  condition  is  that  a  should  be  continuous 

x 

at  the  elastic-plastic  boundary  on  y=0.  By  the  use  of  (4.2a)  and 
(4.3)  we  find 


k-  (ft)  ki  (3  -  ir) 

n  n ' 


(4.9) 


From  (4.8)  and  (4.9)  it  follows  that 


P-.l 

RP  3 

Hence  (4.9)  yields 


(4.10) 


(4.11) 


Equations  (4.8)  and  (4.11)  show  why  we  have  taken  elastic  fields  for 
a  notch  rather  than  for  a  crack.  For  an  elastic  crack-tip  field  the 
conditions  of  reaching  the  yield  condition  at  the  elastic-plastic 
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boundary  would  conflict  with  the  condition  of  continuity  of 
as  can  be  checked  by  setting  p  =  0  in  (4.8)  and  (4.9). 

For  small  values  of  0  and  \fi  we  next  consider  the  continuity  of 

o  and  oD  (i.e.,  the  shear  stress  and  the  radial  stress  in  the  R ,\p 
Riff  R 

system) .  We  use 


o„,  *  (coszil)  -  sin2ii/)o  +  (a  -a  )simJ>cost^ 
R<p  r  xy  y  x 


(4.12) 


and 


o  *  2o  sin^cosip  +  o  cos2ij/  +  a  sin2i p  , 

R  xy  x  y 

in  conjunction  with  (4.2a,b,c)  and  (4.3)-(4.5).  It  may  be 

verified  that  a  is  continuous  to  first  order  in  i|j.  The  stress 
Ry 

a  is  continuous  to  order  unity  by  virtue  of  Eq.(4.11).  Collecting 
K 

terms  to  order  02  and  \p2  yields  the  relation 


(4.13) 


0  -  vl>  ,  where  y  *  /l/2.  (4.14) 

This  result  Implies  that  for  small  values  of  0  and  ^  thfe  center  of 
the  elastic  field  E  is  located  at 

x_  ■  (l-y)x  ,  and  thus  R  ■  yx  (4.15a,b) 

E  p  P  P 

For  y/x  <<  1  it  follows  from  (4.2a,b,c)  that 


-  Ml  -  K>’>  ♦  0<X>* 

(4.16) 

o  -  2k  +  Oft" 

(4.17) 

y  x 

Txy  *  ■*$’  + 

(4.18) 

From  the  stress-strain  relations  (2.4)-(2.6)  we  find  by  using 
(4. 16)- (4. 18) 
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D  ■=  -a(t)|  -(*)*  -  hk(Z)2  +  0(^)4 
x  E  xx  x  x 

D  -  a(c)£!  V)!  +  Ak[l  +  i(*M  +  0(i)* 

y  E  x  x  2  x  x 

D  =  -a(t)  -(^) 3  -  Ak(^)3  +  0(^)5 

xy  E  xx  vx  x 

Now  consider  the  following  expansions  for  y/x  «  1 


u(x,y,t)  =  UQ(x,t)  +  u2(x,t)y2  +  0(y4) 


v(x,y,t)  =  v1(x,t)y  +  0(y3) 


A(x,y,t)  =  AQ(x,t)  +  A2(x,t)y2  +  0(y4) 

The  strain  rates  follow  from  (4. 22) -(4. 23)  as 

3u  3u_ 

D  *»  — -  +  — -  y2 
x  3x  3x  y 


D  =  v, 

y  l 


3v, 


D  »  (2u_  +  "r — ■— )y 
xy  2  3x  7 


(4.19) 

(4.20) 

(4.21) 

(4.22) 

(4.23) 

(4.24) 

(4.25) 

(4.26) 

(4.27) 


By  equating  terms  of  the  same  order  in  (4 ,19)-(4.21)  to  the  corres¬ 
ponding  ones  in  (4 .22)-(4.24)  we  obtain 


3u 
_ _ c 

3x 


0  . 


9u, 

_ t 

3x 


-  k  a 
"  E  ? 


V1  ■  Aok 


3v. 

2“2+  3T 


(4 .28a, b) 
(4.29) 


(4.30) 


By  combining  (4.28),  (4.29)  and  (4.30)  we  obtain 


.  3  v  v.  .  • 

1  _ 1  1  ,  k  a 

2  3X2  7‘  J  E  r 


The  general  solution  to  (4.31)  is 


°l  ■  I  <-2  f  +  ^  +  c<t),«>  . 

p 

and  (4.30)  subsequently  yields 


u2  *  |  [l-^n(f-)l  +  \  ^  -  C(t)x) , 

p 


(4.31) 


(4.32) 


(4.33) 


while  it  follows  from  (4.28a)  that 
uo  -  (k/E)  D(t) 

The  functions  B(t),  C(t)  and  D(t)  must  be  determined  from  continuity 
of  the  particle  velocity  at  the  elastic-plastic  boundary. 

By  the  use  of  (3.10)  we  obtain  from  (4.6)  and  (4.6)  at  the  elastic- 


plastic  boundary 

G«  (k/E)  [U  (t)  +  U,(t)e2]  +  0(0**)  (4.35) 

p  O  L 

Gp-  (k/E)V1(t)6  +  0(93),  (4.36) 

where  (4.10),  (4.11)  and  (4.14)  have  been  used.  In  (4. 35) -(4. 36) 

uo(t)  -  f [ (<y -  +  f>*p "  i(K  ~  (E/w)  (a-3/: 

"  I4  T7  ^  "  2y  -  2<y  +  K+5)xp  +  (>c+5)a]  (E/y)  (4. 30) 
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Vl(t)  "16  y  (lC+1)(V®)(E/W)  (4‘39) 

By  comparing  (4.32)- (4.34)  to  (4.35)-(4.36)  we  conclude 

B(t)  •=  B.a(t)  +  B0x  (t)  (4.40) 

1  2  p 

C(t)  -  [C.a(t)  +  C,i  (t)]/[x  (t)]3,  (4.41) 

1  2  p  p 

where 

B1  "  32  T7  ltc+5  +  AY(<+l)](E/p)  -  |  (4.42) 

B2  =  "32  Y2  ^K+5  +  2y(k+1)](E/u)  (4.43) 

C1  -  32  [  - (k+5)+  2y(k+1)](E/p)  +  |  (4.44) 

C2  "  32  y*  1~(k+5>  +  (k+1) ] (E/y)  <***5) 


In  the  plane  of  the  crack  we  have  ey  ■*  v^.  At  the  elastic  plastic 
boundary  (4.7)  yields  for  small  y 


(eyJPB 


3  k  1,  1. 

8  7  7(K-  3> 


where  (4.11),  (4.14)  and  0  ~  y/x  have  been  used.  In  the  stationary 
coordinate  system,  (4.32)  can  now  be  integrated  to  yield  the  total 
strain  as 


.  .  x  -a(t) 

+  B2  xP/a)*n  -r&r 


.  r  x  -a(t)-|2  k 

r(xi,t)  - 

£  • 

/'  .  .  .  x.-a(t)  x 

{l2  -E  +  B2  ^(xp/a)Hn  -^5-  -  ip/i)  dt 

X  P  P 

P  P 

I  IV 


+l 


.  [x.-a(t)]2 

_•  +  Vp1  K  (t)T  dt 

P  P 


(4.46) 


(4.47) 


Here  t  is  the  time  that  the  elastic-plastic  boundary  arrives  at 
P 

position  x,,  i.e.,  a(t  )  +  x  (t  )  ■  x, . 

1  P  P  P  1 

For  a  given  external  load,  we  presumably  know  in  terms  of  the 

increase  in  "crack  length"  a(t)  +  (l-y)Xp(t).  A  relation  between 

Xp(t)  and  a(t)  can  subsequently  be  obtained  by  the  use  of  (4.15b)  and 

(4.11).  Hence,  in  principle,  a(t)  is  the  only  unknown  quantity  in 

(4.47).  An  equation  for  a(t)  and  a(t)  can  be  obtained  from  (4.47) 

by  the  use  of  the  critical  strain  criterion  for  crack  propagation. 

This  criterion  stipulates  that  crack  growth  will  proceed  when  a  critical 

strain  level  e  is  maintained  for  e  ,  in  the  plane  of  the  crack  at  a 
cr  y 

characteristic  distance  A  ahead  of  the  crack  tip.  It  appears, 
however,  that  it  will  be  very  difficult  to  solve  a(t)  and  a(t)  from 
the  integral  equation  that  can  be  extracted  from  (4.47). 

It  is  of  interest  to  note  that  v^  as  defined  by  (4.32),  (4.40) 
and  (4.41)  contains  terms  proportional  to  a(t)  and  terms  proportional 
to  Xp(t).  The  former  correspond  to  crack  propagation,  while  the 
latter  are  equally  valid  for  a  stationary  crack  under  monotonic 
loading.  By  setting  a(t)  =  0,  we  obtain 


S,1  .  k  ,B2  .  V  , 

sT  •  B  (r  +  TT7E7TT  1  dt  • 


(4.48) 


This  result  shows  that  for  a  stationary  crack  and  plane  stress,  the 


strain  ey  has  a  singularity  of  order  1/x  at  the  crack  tip,  in  agreement 
with  the  results  of  Rice  (1973) . 


17 


The  solution  for  in  the  plane  of  the  crack  simplifies  con¬ 
siderably  for  two  special  cases.  The  first  of  these  is  the  case  of 
a  self-similar  field.  A  self-similar  solution  materializes  when 

a(t)  is  proportional  to  x  (t), 

P 

a(t)  *  Ap  xp(t)  (4.49) 

Substitution  of  (4.49)  into  (4.47)  and  subsequent  introduction  of  the 
new  variable 


s  -  [x  -a(t)]/x  (t)  -  fx  -A  x  (t)]/x  (t)  (4.50) 

1  p  1  P  p  p 

reduces  (4.47)  to 

ey(x/V  "  (V pb  +  EUn(lTTo)1  ~  E(Bi+B2/APHn(orTt)) 

x/x  x/x 

-!{  i2'”")  -  <WV>  Tsr  -!  /  P  <ciAp«2>  Is5-  «-51> 

1  P  1  H  P 

where  we  have  used  that  x  ■  -  a(t).  It  is  noted  that  the  strain 

comes  out  to  depend  only  on  the  self-similar  variable  x/x  (t). 

P 

For  the  second  special  case  all  fields  are  assumed  to  be  time- 

invariant  to  an  observer  traveling  with  the  crack  tip.  This  is  the 

steady-state  case  when  depends  on  x  ”  x^  -  a(t)  only.  Now  we  have 

that  a  ■  i  ■  constant  -  c_,  and  (’)  =  -  c_3  .  Equation  (4.32) 
p  r  r  x 

becomes 


I  <1  -  ;<W  ’ 


(4.52) 
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where  and  C 2  are  defined  by  (4.42)-(4.45) ,  and  ■  constant. 

Equation  (4.52)  may  be  Integrated  to  yield 

%<*>  “  (0PB  +|{Un(~-)]2  -  (B1+B2)S,n(^— ) 

P  P 

■  i(ci'K:2)[(f“)3-11}  (4*53) 

P 

Equation  (4.53)  offers  better  possibilities  for  the  application 

of  a  critical  strain  criterion.  If  e  ■  e  at  x  *  A,  (4.53)  can  be 

y  cr 

solved  for  the  constant  value  of  x^ .  The  corresponding  external  load 
can  subsequently  be  computed  by  the  use  of  (4.15b)  and  (4.11). 
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Fig.  1:  Moving  crack  tip 


Fig.  2:  Center  of  elastl 
plastic  boundary 
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